Topological susceptibility in Lattice QCD with unimproved 

Wilson fermions 



Abhishek Chowdhury, Asit K. De, Sangita De Sarkar, A. Harindranath*, Santanu 

Mondal, Anwesa Sarkar 

Theory Division, Saha Institute of Nuclear Physics 
1/AFBidhan Nagar, Kolkata 700064, India 

Jyotirmoy Maiti 

Department of Physics, Barasat Government College, 
10 KNCRoad, Barasat, Kolkata 700124, India 



Abstract 

We address a long standing problem regarding topology in lattice simulations of QCD 
with unimproved Wilson fermions. Earlier attempt with unimproved Wilson fermions 
at j3 = 5.6 to verify the suppression of topological susceptibility with decreasing quark 
mass (m q ) was unable to unambiguously confirm the suppression. We carry out system- 
atic calculations for two degenerate flavours at two different lattice spacings (/3 = 5.6 
and 5.8). The effects of quark mass, lattice volume and the lattice spacing on the span- 
ning of different topological sectors are presented. We unambiguously demonstrate 
the suppression of the topological susceptibility with decreasing quark mass, expected 
from chiral Ward identity and chiral perturbation theory. 
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1. Introduction 

Because of the explicit violation of chiral symmetry by a dimension five kinetic 
operator, there have been persistent concerns about the Wilson formulation JT] 13 of 
fermions on the lattice in reproducing the chiral properties of continuum QCD. To 
address these concerns, for the past few years, we have been studying (3] O the chiral 
properties of Wilson lattice QCD. We have studied the emergence of the chiral anomaly 
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with unimproved and improved Wilson fermions, the associated cutoff effects and the 
approach to the infinite volume chiral limit in the context of flavour singlet Axial Ward 
Identity to order &(g 2 ). In this work we address various issues associated with topo- 
logical charge (Q) and topological susceptibility (x) in lattice QCD simulations with 
two degenerate flavours of unimproved Wilson fermions. The detailed account of low 
lying spectroscopy and autocorrelation studies will appear separately. 

Earlier attempt [5| with unimproved Wilson fermions and HMC algorithm per- 
formed simulation at j3 = 5.6, lattice volumes 16 3 x 32 and 24 3 x 40 and m K > 500 
MeV. Their results as well as the results from other collaborations [6| were presented, 
adopting a mass-dependent renormalization scheme, by scaling topological suscepti- 
bility and m\ by appropriate powers of quark mass dependent r^ja where ro is the 
Sommer parameter. Since ro/a significantly increases with decreasing quark mass, the 
suppression of topological susceptibility may be concealed in such a plot especially 
for large pion masses and for lattice actions which have more severe cutoff artifacts. 
Results of Ref. [5| failed to show the suppression of topological susceptibility un- 
ambiguously when data was presented in this manner: see Fig. 12 in Ref. J5]. (Note 
however that, at present, a mass-independent renormalization scheme appears to be 
preferred in the lattice literature where ro in the chiral limit is used to scale the data.) 
Since topological susceptibility is a measure of the spanning of different topological 
sectors of QCD vacuum, the inability to reproduce the predicted suppression may raise 
concerns about the simulation algorithm and the particular fermion formulation to span 
the configuration space correctly. Unimproved Wilson fermion has <^(a) lattice artifact 
and hence it is also important to study the effects of scaling violation. We perform 
a systematic study using unimproved Wilson fermions at different volumes, different 
lattice spacings (j3 = 5.6 and 5.8) and m K > 300 MeV. We present our data in both mass- 
dependent and mass-independent renormalization schemes. In the former scheme, our 
data (at both couplings /3=5.6 and 5.8) clearly show the suppression of topological sus- 
ceptibility in the region of pion mass m K < 500 MeV. In the latter scheme, our data 
(at both couplings j3=5.6 and 5.8) exhibit the suppression of topological susceptibil- 
ity in the whole region of pion masses studied. When plotted using mass independent 
scheme, the data of Ref. [5] also clearly show the suppression of topological suscep- 
tibility for the avaiable range of their data (m K > 500 MeV). Since our data exhibits 
suppression in the lower pion mass region independent of the renormalization scheme 
used, we are able to demonstrate unambigously the suppression of topological suscep- 
tibility with decreasing quark mass, expected in continuum QCD in the region of lower 
pion mass. Advancements in both algorithm and technology have made this study 
possible. 

2. Measurements 

We have generated ensembles of gauge configurations by means of HMC Q [8] 
and DDHMC |9| algorithm using unimproved Wilson fermion and gauge actions with 
Hf = 2 mass degenerate quark flavours. At j3 = 5.6 the lattice volumes are 16 3 x 32, 
24 3 x 48 and 32 3 x 64 and the renormalized physical quark mass (calculated using axial 
Ward identity) ranges between 15 to 100 MeV (MS scheme at 2 GeV). At j3 = 5.8 the 
lattice volume is 32 3 x 64 and the renormalized physical quark mass ranges from 20 to 
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Figure 1: Topological charge of gauge configurations versus HYP smearing steps for ]3 = 5.8 and K = 
0.15475 (% ~ 300MeV) at lattice volume 32 3 x 64. 



90 MeV. The lattice spacings determined using Sommer parameter at /3 =5.6 and 5.8 
are 0.077 and 0.061 fm respectively. The configurations for lattice volumes 24 3 x 48 
and 32 3 x 64 and for 16 3 x 32 at K = 0.15775 are generated using DDHMC algorithm. 
The number of thermalized configurations ranges from 2000 to 12000 and the number 
of measured configurations ranges from 70 to 500. 

For topological charge density, we use the lattice approximation developed for 
St/ (2) by DeGrand, Hasenfratz and Kovacs IfTol , modified for SU (3) by Hasenfratz 
and Nieter [11] and implemented in the MILC code [8 1. It uses ten link paths described 
by unit lattice vector displacements in the sequence {x,y,z, —y, —x,t,x, — t , — x, —z] and 
{x,y,z, —x,t, —Z,x, —t, —x, —y} plus rotations and cyclic permutations. To suppress the 
ultraviolet lattice artifacts, smearing of link fields is required. The link field is smeared 
by 20 HYP smearing steps and optimized smearing coefficients a = 0.75, (X2 = 0.6 and 
(Xq = 0.3 lfl2ll . We have observed that smearing brings the topological charges close to 
integer values and behaviour is better for the smaller lattice spacing, as expected. 

In Fig. [T]we show the behaviour of topological charge of gauge configurations with 
smearing steps for j3 = 5.8 and K = 0.15475 at lattice volume 32 3 x 64. It is evident 
that the topological charges of different configurations are clustering about the integer 
values after about 10 smearing steps and values are stable under further smearing steps. 
In Fig. [2] we present the behaviour of topological susceptibility with smearing steps 
for the same lattice parameters which shows that the susceptibility is very stable with 
smearing steps after 10 steps. 

In Fig. [3] we show the Monte Carlo time history of topological charge for /3 = 5.6 
and 5.8 for the smallest and the largest K and there is some evidence of trapping of the 
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Figure 2: Topological susceptibility versus HYP smearing steps for j5 = 5.8 and K = 0.15475 at lattice 
volume 32 3 x 64. 





Figure 3: The Monte Carlo trajectory history for topological charge with unimproved Wilson fermion and 
gauge action for (a) j8 = 5.6 with a gap of 24 trajectories between two consecutive measurements and (b) 
P = 5.8 with a gap of 32 trajectories between two consecutive measurements. 
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topological charge only at j3 = 5.8 and largest K. 

Fig. [4] displays six histograms of topological charge distributions, for two values 
of j3 and different volumes. The topological charge data were put in several bins and 
the bin widths were chosen to be unity centered around the integer values of the topo- 
logical charges for all the cases. From theoretical considerations the distribution of the 
topological charge is expected to be a Gaussian fl3l . Since our configurations are large 
in number but finite, an incomplete spanning of the topological sectors may occur and 
(Q) may not be zero. Hence we define the susceptibility to be 

x = ^((Q 2 )-(Q) 2 )- (i) 

The expected distribution is 

n meas ( Q \ /ON 

nQ= vmm^W) exp vm^W)) (2) 

where n meas is the total number of measurements made. The Gaussian curves in Fig. [4] 
are obtained by using Eq. [2] It is evident from Fig. |4]that for a given j3 and volume 
the width of the distribution decreases as K increases, indicating that the topological 
susceptibility % decreases with decreasing quark mass. 



3. Results 

In table[T]we present, in lattice units, the pion masses, unrenormalized quark masses 
and topological susceptibilities measured. For ease of comparison with earlier presen- 
tations of results and clarification we present for j5 — 5.6 and 5.8 at lattice volume 
32 3 x 64, topological susceptibility versus in the units of ro (quark mass depen- 
dent), (i.e., in mass-dependent renormalization scheme) in Fig. [5] At j8 = 5.6 our 
second largest quark mass (in Fig. [5]) is very close to the lowest quark mass of Ref. 
J5). In this figure the shaded region corresponds to m n > 500 MeV. In the lower pion 
mass region our data for both j3 = 5.6 and 5.8 clearly show the suppression even in 
mass-dependent renormalization scheme. 

Fig. pi shows our results for topological susceptibility versus m\, in the units of 
Sommer parameter (ro) at the chiral limit (i.e., in the mass-independent renormaliza- 
tion scheme), for j3 = 5.6 and at lattice volumes 16 3 x 32, 24 3 x 48, and 32 3 x 64. We 
also show the results of SESAM-T^L collaborations Q. Using the numbers given in 
[5], we have replotted it after scaling by the value of ro quoted at the physical point. 
This new plot clearly shows the suppression of susceptibility with decreasing quark 
mass in the earlier SESAM-T^L data with unimproved Wilson fermion. Our results 
carried out at larger volume and smaller quark masses unambiguously (i.e., indepen- 
dent of the renormalization schemes used) establish the suppression of topological sus- 
ceptibility with decreasing quark mass in accordance with the chiral Ward identity and 
chiral perturbation theory. Further we note that, for a given K, the value of topological 
susceptibility increases with the volume as expected from finite volume considerations. 
This effect is more noticeable in smaller volumes. For large enough volume topological 
susceptibility should be independent of volume, since (Q 2 ) scales with the volume. 
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Figure 4: The topological charge distribution for (a) ji = 5.6, K = 0.15775, volume = 16 3 x 32 (b) /3 = 5.6, 
K = 0.15775, volume = 24 3 x 48 (c) j3 = 5.6, If = 0.15775, volume = 32 3 x 64 (d) J3 = 5.6, If = 0.15825, 
volume = 24 3 x 48 (e) j8 = 5.8, if = 0. 1543, volume = 32 3 x 64 (f) JS = 5.8, K= 0. 15475, volume = 32 3 x 64. 
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Table 1: Pion masses, unrenormalized quark masses and topological susceptibilities in lattice units. 



7 



0.08 



0.06 



0.04 



0.02 



L32T64; p = 5.8 
L32T64; p = 5.6 



■ 



0.5 



1 

2 2 



1.5 



Figure 5: Topological susceptibility versus m| in the units of rg (quark mass dependent) for fi = 5.6 and 5.1 
and at lattice volume 32 3 x 64. 
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Figure 6: Topological susceptibility versus m| in the units of ro (at chiral limit) for j3 = 5.6 and at lattice 
volumes 16 3 x 32, 24 3 x 48, and 32 3 x 64 compared with the results of SESAM-T^L collaborations (5). 
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Figure 7: Topological susceptibility versus m q in the physical units for j3 = 5.6 and j} = 5.8 at lattice volume 
32 3 x 64 compared with the results of mixed action (Clover and Overlap) 1181 . For the latter, the two separate 
lattice spacing determinations, a = 0.0784 fm and a = 0.070 fm quoted have been used. 



Since unimproved Wilson fermion has lattice artifacts it is important to esti- 
mate the scaling violations of our results. Fig. [7] shows topological susceptibility in 
the physical units for j3 = 5.6 and j3 = 5.8 at lattice volume 32 3 x 64 versus nonper- 
turbatively renormalized [14| quark mass in MS scheme lfl5l at 2 GeV. Topological 
susceptibility at j3 = 5.8 approximately matches with that at /3 =5.6 within the error 
bars but the latter data is systematically below the former. This behaviour which is 
qualitatively consistent with leading order lattice artifact lfl6l is observed also by the 
MILC collaboration ifTTl . For comparison, in Fig. [7J we have also shown the results of 
mixed action calculation of Ref . 1 1 8 1 at j3 = 5 . 3 where the authors have quoted two sep- 
arate lattice spacings, a = 0.0784 fm and a = 0.070 fm for the same j3. The calculation 
of Ref. [ 18 1 employs gauge configurations generated with dynamical ^(a) improved 
Wilson fermion, using DDHMC algorithm and topological charge is measured using 
a fermionic operator, namely, the Neuberger-Dirac operator. Fig. [7] shows that our 
results of the topological susceptibility favourably compare with that of Ref. ifTSIl . In 
Ref. |[T9l we have shown that (see Fig. 6 of Ref. |[T9l ) the topological susceptibility at 
/3 =5.8 is consistent with leading order chiral perturbation theory prediction. 

In conclusion, we have addressed a long standing problem regarding topology in 
lattice simulations of QCD with unimproved Wilson fermions. Calculations are pre- 
sented for two degenerate flavours. The effects of quark mass, lattice volume and the 
lattice spacing on the spanning of different topological sectors are presented. By per- 
forming simulations at smaller pion masses we have shown the suppression of the topo- 
logical susceptibility with respect to decreasing quark mass irrespective of the method 
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of presenting the data. The suppression expected from chiral Ward identity and chiral 
perturbation theory is thus demonstrated unambiguously. Furthermore, our results of 
the topological susceptibility favourably compare with that of £?(a) improved Wilson 
fermion. Our results unequivocally show that lattice QCD with naive Wilson fermions 
together with DDHMC algorithm, for the range of quark masses and lattice spacings 
studied, is able to span the configuration space correctly and reproduce the chiral be- 
haviour of continuum QCD. In future work, we would like to explore lighter quark 
masses and smaller lattice spacings. To achieve this goal, however, improvements in 
the algorithm might be necessary. 
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